uE  ET  ?U^ 
/^.     T\  /  H     ^^^  YORK   UNIVERSITY 

^  m  IJl    ^     '"^"*®  °*  Mathematical  Sdences 

.1^      ^^      j>      Division  of  Electromagnetic  Research 
Oc  CCX  X 


RESEARCH    REPORT    No.    EM-70 


Reflection  and  Refraction  of  Electromagnetic  Waves 
by  a  Dielectric  Slab  between  Dielectric  Media 


W.    ELWYN    WILLIAMS 


NEW  YORK  UNIVERSITY 

INSTITUTE  OF  MATHF.M  ATICAL  SCIENCES 

LIBRARY 

25  Waverly  Place,  New  York  3,  N.  Y. 

WAY  5    1955 


CONTRACT  NO.  A  F  -  1  9  ( 1  22)-4  2 
DECEMBER  1954 


NEW  YORK  UNIVERSITY 
Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 

Research  Report  No.   EM-70 


REFLECTION  AND  REFRACTION  OF  ELECTROMAGNETIC  WAVES 
BY  A  DIELECTRIC  SLAB  BETWEEN  DIELECTRIC  MEDIA 

W.   Elwyn  Williams 


W.   ElWyn  Williams 

Morris  Kline 
Project  Director 


The  research  reported  in  this  document  has  been  made  possible 
through  support  and  sponsorship  extended  by  the  Air  Force  Cam- 
bridge Research  Center,   under  Contract  No.   AF-19(122)-42.     It 
is  published  for  technical  information  only  and  does  not  necessarily 
represent  recommendations  or  conclusions  of  the  sponsoring  agency. 

New  York,    1954 


Abstract 

The  general  problem  considered  is  the  reflection  and  refraction  of 
an  arbitrary  electromagnetic  wave  by  a  dielectric  slab  between  two  different 
media.  Explicit  expressions  are  obtained  for  the  radiation  field  in  terms 
of  the  incident  radiation  field.  These  expressions  are  valid  when  the  point 
of  observation  is  not  too  near  the  plane  interface.  The  method  is  ext.ended 
to  cover  the  case  of  one  surface  deviating  slightly  from  a  plane. 
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1,   Introduction 

The  object  of  the  present  note  is  to  present  a  unified  method  of  treating 
the  problem  of  the  reflection  and  refraction  of  an  arbitrary  electromagnetic 
wave  at  interfaces  separating  media  of  different  dielectric  properties.  Several 
special  problems  of  this  type  have  been  considered  in  detail  by  various  authors, 
the  first  complete  solution  being  the  classical  one  of  Sommerfeld  <- J  for  the 
field  of  a  vertical  dipole  above  a  plane  earth.  Horschelmann  >- -> ,  using  the 
same  method  as  Sommerfeld,  obtained  the  solution  for  a  horizontal  dipole  above 
the  earth.  Lurye  •-  J  has  considered  a  more  general  problem  in  calculating  the 
field  of  a  dipole  above  a  dielectric  slab  resting  on  a  plane  earth  of  a  large 
but  finite  conductivity.  A  two-dimensional  version  of  this  latter  problem  had 
been  considered  by  Tai  L  -■  who  made  the  assumption  that  the  conductivity  of  the 
earth  was  infinite. 

All  the  above  problems  may  be  regarded  as  special  cases  of  the  more  general 
problem  of  an  arbitrary  electromagnetic  wave  incident  on  a  dielectric  slab  be- 
tween two  different  dielectrics.  This  general  problem  will  be  considered  in 
the  following  work.  Clearly  it  will  not  be  possible  to  enter  into  a  detailed 
discussion  of  the  properties  of  the  solution  in  a  manner  comparable  to  that  of 
the  above  authors,  and  no  attempt  is  made  to  do  this.  It  is,  however,  possible 
to  obtain  some  fairly  simple  results  concerning  the  form  of  the  solution  at 
large  distances  from  the  source. 

The  method  used  is  an  extension  of  that  adopted  by  Copson'-  J  for  considering 
the  diffraction  of  a  plane  harmonic  electromagnetic  wave  by  an  aperture  in  a  per- 
fectly reflecting  plane.  The  components  of  the  electric  field  wMch  are  in  a 
direction  tangential  to  the  plane  interfaces  are  constructed  so  that  they  are 
continuous  across  both  plsmes.  The  components  of  the  magnetic  field  are  then 
expressed,  using  Green's  functions,  in  terras  of  the  tangential  electric  compon- 
ents on  the  interfaces.  Integral  equations  are  then  obtained  for  these  latter 
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components  by  applying  the  condition  of  continuity  of  the  tangential  components 
of  the  magnetic  field  across  the  planes.  These  integral  equations  are  then 
solved  by  means  of  double  Fourier  transforms. 

The  method  is  also  extended  to  cover  the  case  of  one  surface  deviating 
slightly  from  a  plane.  The  only  difference  in  this  problem  is  that  the  boundary 
conditions  of  continuity  are  now  to  hold  on  a  surface  which  is  not  quite  plane. 
The  method  adopted  is  tc  consider  the  boundary  conditions  on  the  plane  from 
which  the  actual  surface  deviates  slightly.  The  effect  of  this  approach  is 
that  the  components  tangential  to  this  plane  will  now  have  known  discontinuities 
across  the  plane*,  the  formulation  will  otherwise  be  the  same  as  for  plane  surfaces, 

The  general  results  for  a  nearly  plane  surface  will,  of  necessity,  be  quite 
complicated,  but  simplification  is  afforded  by  considering  speciad  fonns  of  the 
surface  and  of  the  incident  field.  In  any  particular  problem  it  should  be 
possible  to  adopt  the  method  used  by  Rice  *-  -^  in  considering  a  perfectly  reflecting 
nearly  plane  surface  and  compute  various  averages  of  the  field  components,  the 
average  being  taken  over  many  different  surfaces  possessing  the  same  statistical 
properties. 

2.   Mathematical  formulation 

It  will  be  assiimed  that  there  are  three  dielectric  media  and  that  they 
occupy  the  regions  1,  z  >  Oj  2,  -£.  <  a  <  Oj  and  3,  z  <  -a,  respectively, 
where  z  denotes  the  normal  to  the  planes  (see  Fig.  1).  The  x  and  y  axes  of 
a  system  of  Cartesian  coordinates  are  chosen  parallel  to  the  planes  in  such 
a  manner  that  0(x,y,z)  is  a  right-handed  system  of  axes.  For  simplicity 
it  is  assumed  that  the  magnetic  permeability  is  the  same  for  all  the  media, 
and  a  time  factor  e    will  be  understood  throughout  the  work. 
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Figure  1 


Maxwell's  equations  in  each  of  the  three  regions   (in  normalized  M.K.S. 
units)  are 


(1) 


curl  H  =   ( (T"    +  iooe.e  )E,  carl  E  «  -icd(j,  H,         J   »  1,2,3, 


where  H  is  the  magnetic  field,  F  the  electric  field,    (y.  the  conductivity,   ard 

6.  the  dielectric  constant;   6   ,  u.     are  the  values  of  the  dielectric  constant 
j  *     o*  ^o 

and  magnetic  permeability  in  vacuo.     It  follows  from  equation  (1)  that  all 
the  Cartesian  components  of  E   and  IT  satisfy  Helmholtz's  equation 


(2) 


(V^  +  k^)  Y  -  0, 
J 


J  =  1,2,3, 


2  2 

where  V     is  the  Laplacian  operator  and  k.  ■  -ifon  (  (p.+ioac  .e   ) ,     Throughout  the 

J  J  J 

work  it  is  assumed  that  k.  has  a  negative  imaginary  part  (i.e.,  that  the  con- 
ductivity  never  vanishes) j  this  assumption  gives  a  simple  form  for  the  radia- 
tion condition  at  infinity,  namely  that  all  the  components  vanish  exponentially 
with  large  distances.  It  also  ensures  that  all  th3  integrals  occurring  in  the 
analysis  are  convergent. 
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On  applying  the  method  adopted  by  Copson  it  is  fairly  simple  to  obtain 
the  following  integral  expressions  for  the  Cartesian  components  of  E  and  H  in 
the  region  z  >  0,  The  subscripts  x,y,z  in  the  formulas  denote  the  component 
direction- 

(3) 


3U-, 
^x--3z^       ^y^ 

■  az    » 

au-j^     av^ 
^z  '  ax   *  ay   ' 

^- 

^1^1     ^1 

\' 

where 

"i "  h/^1  ^x'^'^y*  ^1 "  ^  Ai  ^ydx'<^y'»  \  '  hi/h  ^z'^'^y- 

Here  (jS  '  e'    l/r,  where  r     =  (x-x')     +   (y-y')     +  z   >  and  the  script  letters 
denote  the  values  of  the  respective  components  on  the  plane  z'    =0.     Unless 
otherwise   stated  all  integrals  will  be  xinderstood  to  be  taken  over  the  whole 
(x'y '  )-plane.     By  Maxwell's  equations,  the  quantities  U-.,   V^,  W^  are  related  by 

1  1^1      \    ax  dy 

Equations  (3)  apply  for  the  region  z  <  -a  with  only  the  difference  that  the 
signs  of  the  expressions  must  be  reversed  and  k,  replaced  by  k^.  This  gives 


for  z  <  -a 


(U) 


au^                   av^                      au,  av^ 

X     az    »          y      az    *        ^z         ax  ay    » 

k^v^     aw                 k^v^  3W^  aw 

H      -    -   T^-i  +   ^    ,         H      -   r-^-^  +  ^   ,  H„    -   ^    , 

X          iMti         ax    *         y       icou.  ay    '  z       az     ' 
^0                          «                   o 


where 


"3  "  2^  A3^x^^'^y'»      ^3  "  7W  A3^/^'^y'      ^3  "  7^ /'^3'^z'^"^' ' 
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Here  the  primes  denote  values  on  z'  -  -a  -  0,  and  ^^   «  exp(-ik-r  /r  ),  where 

2        2       2      2 
r»   -  (x-x')  +  (y-y)  +  (z+a)  . 

In  order  to  calculate  the  field  components  in  the  region  -a  <  z  <  0  it  is 

necessary  to  know  the  Green's  functions  G,(x-x«,  y-y*,   z,z'),  G^Cx-x' ,  y-y',   z,z') 

for  the  region  -a  <  z  <  0,    |x|,    |y|  <   oo.     These  are  defined  as  solutions  of 

(V^  +  k2)G  -  -Un6(x-x')6(y-y')5(z-2') 

such  that  G^  -  0  for  z'  =  0,  -a  and  BGp/Bz  =  0  for  z'  ■  0,  -a;  5(x-x")  is  the 
Dirac  delta-function.  Both  these  functions  have  the  common  property  of  behaving 
like  diverging  waves  at  a  large  distance  from  (x'y'z').  The  exact  determination 
of  the  Green's  functions  is  considered  in  the  following  section. 

It  is  easy  to  show  that  the  field  components  for  region  2  are  given  by 

(5)       E.--r7/^^;'feL^-  ^x  (^)z.--[^-'<^y'> 


"/*^l\       ^m/^ii 


(6)      K  '  '  n^/  CU^  z'=0  -  ^y  I  3^  z'.-a   ^'^^^ 


(^)      H^-Hn/   -t4  h)z.^^'-(^2)z.«-a^r][^-'^' 


k!    /> 


(8)      H  -  ^/  [(G,)    ^;-(G,)    c]^-'^' 

•'     ^o/  *—   z'"0      z'  =  -a  -^ 


» 


where  the  triple  primes  denote  values  on  z  =  -a+0  and  the  double  primes 
denote  values  on  z  =  0  -  0.  The  expressions  obtained  for  H^  and  Hy  by  appli- 
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cation  of  Green's  theorem  to  G^   and    /i/     and      /     will  actually  involve 
9  ^^92,  3/7_/3i>     but  substitution  from  Maxwell's  equations  for  these  ternis 
and  some  algebraic  manipulation  will  yield  equations   (7)   and  (8),     It  should 
also  be  noted  that  the   A^   s  in  these  equations  may  be  expressed  in  terms  of 
X',  y' -derivatives  of  the  tangential  electric  components. 

If  there  is  an  incident  field  E  ,  H     in  the  region  z  >  0,  then  the  total 

i         r         i        -T  t       t 

field  may  be  written  asE     +E,     H     +lrinz>0  and  as  E  ,  H     in  z  <  0. 

It  is  necessary  that  F     and  E     be  continuous  across  z  -  0,   and  hence 

X     y  * 

If  we  write  t      and  £   for  E  ,  E  on  z  =  0  then  E  ,  E^  at  any  point  will  be 
X      y     x'y  x'y     ^   ^ 

given  by  equation  (3).  The  part  of  the  integral  in  equation  (3)  involving 

-E  ,  -E^  on  z'  =0  clearly  corresponds  to  the  reflected  field  in  the  region 
X    y 

z  >  0  when  the  plane  is  perfectly  reflecting.  Hence,  if  we  use  the  super- 
script 0  to  denote  the  field  in  z  >  0  when  the  plane  is  perfectly  reflecting, 
/='  "  „   „t  „t 

X  '  ^ 

are  given  by 


and  write  ^  "  ,   ^        for  E  ,  E  on  the  plane,  the  field  components  for  z  >  0 
X    y      X   y 


(9) 
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•         ^if     ^1         „iti 


If  vre  put  ^      '^  S   >   £ .,  "    ^  „*  then  equations   (U)   -  (9)  represent  the 
X  X         y  y 

solution  of  Maxwell's  equations  in  all  three  regions  with  E  ,   E     continuous 

X   y 

across  z  =  0,  z  »  -a.  It  is  clear  that  the  condition  of  continuity  of  H  ,  H  , 
across  the  surfaces  will  give  rise  to  four  integral  equations,  but  the  actual 
equations  will  not  be  given  here.  In  these  equations  the  terms  (H  ),  (H  )  on 
z  =  0  will  occur  and  these  two  terms  must  be  replaced  by  2H  ,  2H  on  z  =  0. 
The  solution  of  these  equations  is  considered  in  Section  U. 


3.   The  determination  of  the  Green's  functions 
Both  G, ,  Gp  are  solutions  of 

(V^  +  k2)G  -  -lin5(x-x')5(y-y')6(z-z'). 
The  double  Fourier  transform  g  of  G  is  defined  by 

00 

(IC)  S  =  n  /  y  G(t^,t2,z,z')exp[-i(at^  +  pt^^dt^dt^  . 

2        2 
The  behavior  of  G  as  (x-x')  +  (y-y')  ->oo  ensures  that  the  integral  of 

equation  (10)  converges  for  all  real  c  and  p.  We  also  know  that  the  differ- 
ential equation  for  g  is 

i-|  -  ^l^  C  -  -2  5(2-z') 
dz 

where  y,-  =  a     +  p     ~  ^o>   ^^^  S  ^®  continuous  at   z  =  z',  but   (-r^)  -   (t^.)  ^ 

c.  c  dz  ,   t^       dz  ,   ^ 

z='Z'-0  z=z'+0 

Thus  we  have  the  following  expressions  for  g, ,   gp! 

p     sinh  HpZSinh  tip(z'+a)  „     cosh  ix_(z'  +  a)cosh  ^pZ 

^1  tip  sinh  tipa  *  ^2     ^p  sinh  ^p  a  *       -       * 

P     sinh  tipZ'sinli  iip(z+a)  P     cosh  p,p(z+a)cosh  tip^' 

^1  "  ~  jir  sinh  |Xp  a  *  ^2"  ilT  sirli  tip   a  *       — 
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¥e  take  that  branch  of  the  root  [i.   which  reduces  to  ikp  for  a  »  p  =  0»  It  is 
also  known  that  the  double  Foxirier  transform  of  (jf,   »  J^  may  be  written  as 


ll"  ^  exp|-|z-2'|n-Lj  . 


U.   Solution  of  the  integral  equations 

All  the  integral  equations  occurring  in  our  problem  are  of  the  form 

(11)  f(x,y)  -  ^  /g(x',y«)h(x-x',y-y«)dx'dy«. 

Denoting  the  double  Fourier  transforms  of  the  above  functions  by  capital  letters, 
we  may  write  equation  (11)  as 

F  -  GH, 

as  is  well  known.     Now  treating  the  integral  equations  occurring  in  the  present 
analysis  in  this  manner,  we  obr_,ain  the  following  equations: 

(12)  2ia^^H*^  =  g[(a2-kj)  i^  *  J^  ('=^^-4^"^^^(  li  *  h^  *  §1  lu^^^-a^) 

+  tip  iu  ^1  » 


(13)  ^i^gHy"-  =  f 


J^(k^.p2)  ^  l^^^rP^)]  *  «Pg(  ll  *  IP*  ^1  U^^'K 


-  aP  iu  gl   > 


(m         gijOc^-a^)  I2  -   (k^-a^)  13]   *  apfi[j2  *  I3]  -  ^^l'^^  Ik  *  °^^  Ih 


(15)    f. 


i[(p^-4)  h  *  (P^-4^  53]  -  '^^^i[52 "  13]  ■  ^  iil^^'4^  -  °^  h 


where  J^   =  l/ti^,  ^^  -  l/n^,  J^   =  (l/n2)coth  ti^a,  |^  -  (l/n2)cosech  ^ga,  and 
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H*^,  Hy"-,  f,  g,  f^,  g^  are  the  Fourier  transforms  of  (H^)^,  (Hy)^,  f  ^',  ^^,  Ty  . 

The  solution  of  equations  (12)-(1^)  has  been  obtained  explicitly,  but  the 
rather  complicated  expressions  obtained  are  not  given  here.  In  a  later  section 
it  is  shown  that  the  far  field  depends  on  the  values  of  the  functions  at  one 
point  in  the  op-plane,  and  it  is  possible  that  the  best  method  of  obtaining 
this  far  field  is  to  solve  equations  (12)-(15)  numerically  for  the  particular 
values  of  a  and  p  which  are   le quired. 


5.   The  solution  of  special  problems 

In  this  section  the  results  of  the  previous  sections  will  be  applied  to  the 
problems  considered  by  Sommerfeld  and  Horschelmann. 

(a)  Sommerfeld ' s  problem 

In  this  problem  a  and  k-  approach  oo.  The  incident  field  is  due  to  a 
vertically  polai'ized  Hertzian  oscillator,  whose  Fourier  transforms  are 
easily  found  to  be     „  2 

.-i    h^         -h.^      „*i      ^       -^^ 

vhere   h  is  the  height  of  the  dipole  above  the  earth,  (z  "0)0 

If  these  expressions  are  substituted  in  equations  (12)  and  (13)  then 
the  following  expression  is  obtained  for  f : 

2 
f  ■  — 2 ? ® 

Application  of  the  Fourier     inversion     theorem  to  the  function  U-j^  given 
in  equation  (3)  yields 

/2 ok,  tju exp  [-|j- ( z+h )^ exp  [i ( ox+py )3 dadp 
i-t J: . 
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Now  we  make  the  following  substitutions  in  equation  (16):  a  ■  Xcos  6, 
P  «  Xsin9,  y  ■  psin0,  x  ■  pcos^.  Then  on  application  of  the  identity 

cos  <^   J^(Xp)  .  ^  /  e^'^P'^^^^^-^^cos&de 

it  is  seen  that  U-   may  be  rewritten  as 

yr-kf  exp[-  vi^-k?(z+h)lj  (Xp)\dX 
(17)  "         "     '"'"    '  ^  1  J  o 


/?Ik^  (k^  y>?^  +  k^  V^^  ) 


Equation  (17)  agrees  with  the  result  given  by  Sommerfeld. 

(b)     Horschelmann's  problem 

This  problem  differs  from  that  of  Sommerfeld  treated  above  only  in  that 
the  polarization  of  the  Hertzian  oscilator  is  now  horizontal.     We  have 


"x    "  °'      "y    "  ^exp[-M^h]ico^^ 


If  the  expression  for  H   is  substituted  in  equations  (12),  (13)  and  we  pro- 
ceed  as  in  the  previous  problem,  then  an  application  of  the  integral  repre- 
sentations for  J  (Xp),  J2(Xp)  gives 


/?  3  2 
/  X^k^  exp  [-M-  ( z  +h  )j| 

(18)        V"/  — h — T 


jQ(Xp)-J2(Xp)cos2g( 


00 

0      /  J  (V>)exp|^lju(z+h)3 

+  2k;  /  -A h5^ XdX. 

J-  /    \\^-^  V^)\^ 

"^  o 

In  equation  (18)  it  is  assumed  that  p- ,  Hp  refer  to  yx  -k, ,  /X  -k_  . 

Horschelmann  in  his  solution  used  a  Hertz  vector  "ff  of  two  components 
Jf  ,  Jf  ,  and  the  parts  of  fT  reflected  because  of  the  imperfect  conductivity 
of  the  earth  were  fovind  to  give 
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00 

/  J  (Xp)expr-n-  (z+h)] 
(19)    H    '2        -2. !_i ±  XdX, 

op  2 

/  J^(Xp)exp[-tJu(2+hr|lt, 
div  7[=  -2cos^  /  -=—^ ^ i  XdX. 

2  T-r-      3 

It  is  now  easily  verified  that  the  field  given  by  E  =  k.,  T|~  +  ^-  (div  "]]") 
and  derived  from  equation  (19)  agrees  with  that  given  by  -3U, /8z,  where  U-,  is  de- 
fined in  equation  (18) o 

In  this  particular  problem  the  present  approach  hias  an  advantage  over  that 
of  Horschelmann  since  in  his  solution,  although  only  one  component  of  the  Hertz 
vector  was  present  in  the  incident  field,  the  total  field  requires  two  components 
of  the  Hertz  vector.  The  boundary  conditions  at  the  interface  would  otherwise 
lead  to  a  contradiction.  The  present  formulation  avoids  all  these  difficulties. 


6«   The  asymptotic  forms  of  the  field 

In  this  section  we  shall  consider  the  form  of  the  field  at  large  distances 
from  the  slab  (for  this  work  it  will  be  assiomed  that  the  sources  of  the 
incident  field  are  at  a  finite  distance  from  the  plane  z  ■  0)« 

For  convenience  we  shall  replace  x,  y,  z  by  x, ,  x^,  x^.  Hence  we  have  to 

2    2    2    2 

consider  the  form,  as  R  »  x-  +  x_  +  x,  -o  0,  of  the  following  type  of  integral; 

(20)  i-  -^-—   /  9f,  f(xSy')dx'dy  ' 

ax^ax^-  ^ 

where  r  +  /-  »  s . 

For  X',  y'  small. 
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^1 

/v 

e 
R 

exp 

and  also 

3^ 

8x-ax^^ 

I 

-ik^R 

e 
R 

exp 

-ikj^R 

ik^(x^x'  +  x^yO/R  +  0(r"^) 


fik^Cx^x'+x^yO/RJ 


-ik,  R 
ik-  s     .     1      I — 

(-^)  xjxj^^-^^ exp  ik^(x^x'+y^y«)/R 


+  0(R"'). 


2   2 
If  we  assume  that  the  main  contribution  to  (20)  for  R  large  is  from  (x«  +y'  ) 

small,  then  the  asymptotic  form  of  (20)  as  R  ->oo  is 


(ik.)«  x^f  -i'^l^   /   p 
(21)  = —  -i— *!■  ^-^5 /  exp  ik-(x«cos6cos^  +  y'sinOcos^) 

2  TT       R^      -/         ^ 


dx • dy ' , 


where  x_=Rcos©cos^,  x,  =  R  sin9co5^,  x.  =  R  sin(^.  The  integral  in  equation  (21) 
is  seen  to  represent  2n  times  the  Fourier  transform  of  f  at  a  ■  -k^cos6cos0, 
P  •»  -k_sin6cos^. 

If  the  above  approach  is  applied  to  the  following  identity: 


then  we  have 


1   /  ^1  -i 

7^         -§i^H^(x',y«,0)dx'dy'  =  H^(x,y,-2),    z  >  0, 


-ik  R 
ik-sin0  — = H^  (-k-cos9cos(^,  -k,sin9cos^)  =  h  , 


where  h  is  the  radiation  field  of  H  at  the  image  of  the  observation  point  in 
the  plane  2=0.  If  the  incident  field  is  due  to  a  source  distribution  contained 
in  a  finite  volume  then  h  may  be  regarded  as  the  field  at  the  observation  point 
due  to  the  image  of  the  source  distribution  in  the  plane  z  =  0.  Thus  it  is  pos- 
sible tc  express  the  total  radiation  field  directly  in  terms  of  that  of  the  in- 
cident field. 
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The  present  method  of  expanding  the  integrand  gives  the  same  result  as  that  ob- 
tained by  applying  the  saddle-point  method  for  double  integrals  •-' J  to  the  integrals 
of  the  Fourier  inversion  formula. 

Finally  it  should  be  noted  that  the  present  method  is  not  applicable  if 

z  is  small,  i.e.,  if  R  is  large  and  {2  small,  since  in  tnis  case  the  exponential 

in  the  integrand  of  (20)  becomes  small  and  it  is  not  possible  to  neglect  the 

contributions  to  the  integral  from  the  region  near  (x',y')  =  (x,y)« 

For  the  case  of  a  plane  earth  the  radiation  field  due  to  the  imperfect 

conductivity  is 


2         2 
XE  =  -2iocn,  h  k,sin©cos9cos  ^ 


k,sin^ 


,/r75 — 57 


+2ici:n  h 
o  y 


2    2        2  2   2    2—1 
k, (cos  Ocos  ^-1)   (kp-k,co8  ©cos  ^) 

+  


sin  (2 


/kp-k,cos  52 


2         2 

XE  »  -21(41  h  k- sin6cos©cos  Of 
y     ^^o  y  1 


k-sinp 


yl5^ 


cos  (jf 


-2ic4i  h 

^^O  X 


n2  2         2  222         2n 

k^(l-sin  6cos  ?()        (k_-k,Ein  9cos  (2 

—  -«■  ■ 


sin^ 


/7~2 T" 

ykp-kTcos  ^ 


-XH  sin^l        (k?-k^)h  sin9cos©co6^9            |                     fk^sin  ©cos^^i-k^cos  ©cos^0] 
-,5-iL- -^     ^     7 *  h{k.  sin^  -  i-i ••  


-22         222         2^ 
[k,  sin  ©cos  (ji-)>i„cos  ©cos  ^J 


-XH  sln^l 

^5— 2- -  h     <  k-  sinQ!  -  _____ 

°^  '  yk^cos^<2 


-  h 


2     2  2 

(k--k2)8in&cos*cos  (jl 


yk|Ik|Tos^ 


-  lU  - 
where 


(k-sin^  +  ylc|-k^cos^^2  )(  yic|-k^c06^^  k^  +  k2k^sin{2  ) 
sin  9    ykp-k-  cos  ^2 


7 .   Solution  for  one  surface  deviating  slightly  from  a  plane 

In  this  section  we  consider  the  problem  in  which  one  interface,  namely  the 

one  at  z  ■  0,  is  not  plane  but  defined  by  z  =  £h(x,y),  where  e  is  small. 

The  solution  for  a  plane  interface  will  be  denoted  by  the  superscript  (l), 

and  the  total  solution  is  assumed  to  have  the  form 

The  boundary  conditions  to  be  imposed  on  z  =  eh(x,y)  are  that  the  components 
of  the  electric  and  magnetic  fields  tangential  to  the  surface  be  continuous  across 
the  surface.  Hence  to  the  first  order  in  e, 

a  similar  equation  holds  for  E  . 

On  expanding  equation  (22)  in  powers  of  e  and  neglecting  the  terms  of  order 
greater  than  the  first,  we  obtain 

(2)  ^  (2)   3h  ,  (1)    (1).   ah  .  (1)    (1). 

(2) 
and  similarly  for  E   , 

Let  us  denote  by  S^''^^\    £^^'^   the  values  of  E^^  ,  E^^^  on  z  -  0  -  0} 

X      y  X    y 

then  equations  (5;  -  (8)  hold  for  E    etc.  in  the  region  -a  <  z  <  0  if  / 

and  5   are  replaced  by  i  ^  '   and  >:  ^   ,  The  field  components  E  (z)  etc,  in 

7  <   -a,  will  be  of  the  sarne  form  as  before  except  with  the  superscript  (2)  in  each 
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expression.  The  field  components  in  z  >  0  will  now  be  given  by  (3)  with  C 
and   (T  replaced  by 

X         3z  ^  x2     xl  '   ^  zl     z2  ^  az 

and  the  corresponding  expression  for    6    .     The  boundary  conditions  on  H       ,  IT. 
on  z  -  0  are 

(23)  4f  .  h(^'  .  h  ^  (h(,"  -  „ii), 

Equations  (23),  (2U)  contain  no  terms  in  3h/3x,  3h/Sy,  since  by  our  assumption 
of  unit  permeebility  H    is  continuous  at  z  »  0, 

If  the  above  changes  are  carried  out  im  the  formulas  and  the  boundary  con- 
ditions (23),  (2U)  are  applied,  then  a  set  of  integral  equations  will  be  obtained 

(2) 
for  (T    similar  to  those  obtained  in  the  plane  case.  The  only  difference  is 

that  the  incident  field  in  the  latter  equations  will  have  to  be  replaced  by 

functions  involving  the  field  components  for  a  plane  surface.  The  explicit 

form  for  these  functions  as  Faltung  integrals  has  been  obtained,  but  in  view 

of  their  complexity  they  will  not  be  given  here.  Some  simplification  would  be 

afforded  if  the  surface  were  assumed  to  have  a  special  form,  and  it  would  be 

possible,  with  a  considerable  amount  of  labor,  to  carry  out  an  analysis  for 

any  particular  given  field  using  an  approach  similar  to  that  of  Rice, 
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